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ABSTRACT 

We use dimensional reduction techniques to relate real time finite T correlation func- 
tions in (2+1) dimensional QCD to bound state parameters in a generalized 't Hooft 
model with an infinite number of heavy quark and adjoint scalar fields. While static 
susceptibilities and correlation functions of the DeTar type can be calculated using only 
the light (static) gluonic modes, the dynamical correlators require the inclusion of the 
heavy modes. In particular we demonstrate that the leading T perturbative result can be 
understood in terms of the bound states of the 2d model and that consistency requires 
bound state trajectories composed of both quarks and adjoint scalars. We also propose 
a non-perturbative expression for the dynamical DeTar correlators at small spatial mo- 
menta. 
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1. Introduction 



About 10 years ago, DeTar IQ conjectured that the finite T excitation spectrum of 
QCD was much richer than previously expected - excitations with wavelengths longer than 
the magnetic scale l/g'^T are color singlets, at the electric scale 1/gT there are colored 
plasmons and at the hard scale 1/T colored quarks and gluons. He also proposed to test 
this conjecture by studying color singlet current-current correlators of the type 

^""(^) = ^(//^^^(^'^)rX^,^) /^''cir'^(0,r')rXO,r')^ , (1.1) 

where (■ ■ ■) is the average over the Gibbs ensemble. Lattice measurements of these "De- 
Tar correlators" in Euclidean space, relevant for studying screening effects in the finite 
temperature theory, suggest that the screening masses are about 27rT for the light vector 
and pseudovector excitations, and smaller for the scalar and pseudoscalar. 

Some years ago, it was suggested 0, Q that the lattice results on the screening lengths 
and the spatial extension of the static correlators at high temperature could be under- 
stood in terms of a dimensionally reduced theory which is known to confine in the high 
temperature limit. It was shown that the static correlation functions can be calculated 
from an effective theory involving nonrelativistic quarks with colored Coulomb interac- 
tions. The screening lengths were found to asymptote the free quark values, with spatially 
correlated quarks in the transverse directions at all temperatures. These ideas have been 
further developed in several papers, notably those of Huang and Lissia [§, |^, who have 
studied static properties at high temperature in QCD and in the Gross-Neveu model. 

In spite of our increased understanding of the static aspects of high temperature 
QCD, DeTar's main point, namely that the high temperature phase of QCD is confined 
in the sense that all low energy excitations are color singlets, is still neither proven nor 
disproven. The main difficulty is that although there are many sophisticated and efficient 
techniques for calculating Euclidean correlation functions, there are no good methods, 
neither analytical nor numerical, to obtain the real time, finite T correlators. In principle 
these can be obtained by analytically continuing the Euclidean functions, as is done 
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routinely in perturbation theory, but on the lattice this is very difficult without having 
a rather detailed a priori knowledge about the analytic structure of the finite T spectral 
functions (for some alternative attempts in QCD see and in simple model systems see 

There are several lines of attack on this problem. The most direct is perturbation 
theory, where there are well developed methods for calculating real time correlators |]^, 
but since QCD is known to be nonperturbative at large distances even at infinite T, this 
method is of limited value. Another way is to make some specific assumption about 
the spectral functions, as for instance pole dominance with T dependent masses, and 
then use sum-rule techniques to check the consistency of the assumptions and estimate 
spectral parameters in terms of condensate values. This method is likely to be good at 
temperatures well below Tc @, |TI], [l^, but has also been attempted at large T |T^. 



The last approach is to attempt exact, or at least controlled, analytical calculations in 
simple models, with the hope that some of the features can be generalized to realistic 
cases. The present paper is devoted to a study of this kind using 2+1 dimensional QCD 
as a model theory. 

The main idea of our calculation is to express the real time finite T correlators in 2+1 
d QCD in terms of spectral parameters (bound state masses, form factors) of a generalized 
't Hooft model, which is obtained simply by rewriting the theory in 2d language. This 
is of course only advantageous if the heavy {i.e. nonstatic) modes can be totally, or 
to a large extent, ignored in the high T limit. This is known to be the case for static 
observables, and we shall argue that for zero two-momentum, (f = 0, this is true also for 
certain dynamical correlation functions in the sense that only the heavy modes of the 
quark fields have to be retained. This claim is based on a comparison between the high 
temperature limit of perturbation theory in the 2+1 d model and the generalized 't Hooft 
model. This comparison also shows that for nonzero q, the heavy scalars (originating from 
the non-static gluon modes) must be retained, which implies that the two-dimensional 
method is of limited value. However, in this case we can show that consistency requires 



the existence of linear trajectories of bound states containing heavy adjoint scalars in 
addition to quarks and antiquarks. We shall briefly discuss this somewhat surprising 
result in relation to the recent work on QCD2 with adjoint scalars. 

In the next section we show how to reformulate 2+1 d QCD in a two dimensional lan- 
guage, recall some basic results about the 't Hooft model and generalize certain formulas 
to the case of unequal mass bound states. Section 3 is devoted to a study of the x and t 
dependent quark susceptibility which is calculated both in perturbation theory and from 
the 2d bound states. In section 4 we suggest a non-perturbative expression for non-static 
DeTar correlators at small q, and we close in section 5 with a discussion of our results 
and some general comments. 

2. Dimensional reduction 

Our starting point is the Euclidean action for massless QCD3, 

S3E = J d'xCsE{A^J,i) , (2.1) 

where 

C3e{A^, 4^, = ^TtF^^F^p + k-^lJa + 93laAo,)^ , (2.2) 

with the conventions x"' = (x*,r), a,f3 = 1,2,3 and the gamma matrices are defined by 
(7^,7^,7'^) = (— cr^, cr^, cT'^) with cr* the usual Pauh matrices, is the (dimensionful) cou- 
pling constant in three dimensions. Next we simply rewrite the action in two-dimensional 
language, using dimensional reduction techniques. We differ however from e.g. ref. 
in that we keep all the nonstatic modes as well. These heavy modes will then manifest 
themselves as massive fields in the 2d action. We start by making the following mode 
expansion of the quark and gluon fields 

°° . f 

n=—oo 
00 

i"(x\r) = E i^(xOe-^<- , 



where = {2n + 1) vrT and = 2rmTT. Using the static gauge condition r) = 

0, which imphes that the only remaining mode of the field is a (perturbatively) massless 
Higgs A^(a;*, r) = A^{x'^) = we get the following final 2d form of the action, 

+ /3 5:|rf2x^„^737,ir"^^™ + ^2'f AU,"] , (2.3) 



where j = 1,2 and /5 = . 5*1^^^ [A*, 0] is the action for a two-dimensional Euclidean 
YM-Higgs model, as discussed in detail in e.g. i, [H, and S^ilA^.A"]] is the action for 
the nonstatic gluonic modes including the coupling to the static ones. 

As discussed in 0, static DeTar correlation functions can, up to 0(1/T) corrections, 
be calculated from an effective dimensionally reduced theory, although this is not possible 
in 3+1 d QCD. The only effect of the heavy scalar fields is then to renormalize the coupling 
in the 2d model and to add new vertices. One must of course keep the heavy quark fields 
since there are no massless fermionic modes due to the anti-periodic boundary conditions. 
Even this model is too difficult to treat analytically, and we will make the simplifying 
assumption that the scalar field can be neglected altogether. Arguments for why this 
does not change the basic properties of the model were given in ll^ , where it was argued 



that the net effect of the scalar field is to renormalize the string tension and to introduce 
an effective XjR potential at large distances. A strong indirect evidence for the validity 
of this approximation was given in ref. |T^, where it was shown that the leading high 
temperature behavior of the susceptibility is correctly reproduced by the contributions 
from 2d QCD bound states solely, neglecting any effects of the Higgs field. 

When it comes to nonstatic correlation functions the situation is much more compli- 
cated. There is no a priori reason for neglecting the heavy gluon modes, and as we shall 
see below in the case of the non-static susceptibility, in general they do contribute. 

To apply the 2d QCD technology, we follow the approach in [0 and rotate to a fictitious 
2d Minkowski space by: 

t = x^ = -ix^ 82 = -ido Aq = iA2 



with /i, z/ = 0, 1 and the Minkowski metric (H — ). We also perform a chiral rotation on 
the quark fields and a rescaling (in order to get the mass term on the conventional form 
and to give the quark fields the right canonical dimension p[) 

7/>„ = VTe^^'^'V^, (2.5) 



With this, we get the following final form for the action we shall use in the static sector, 

SsE = -iS^f^m-tJ2 fd^xii^i^Y^^-grA,-Mn)i;r^ • (2.6) 

In ( p.6|) the gluon field has been rescaled as well and g = VTg^ is the 2-dimensional 
coupling constant. The model described by (|2.6| ) is QCD in two dimensions with an 
infinite number of fermion fields with masses M„ = {2n + l)vrT, which can be solved 
in the large Nc limit, as originally shown by 't Hooft and subsequently elaborated 



by many others. We will show that the methods developed in e.g. pO| and [0] can 
be adapted to the present problem, and that the difficulties associated with the infinite 
number of quark fields can be handled. 

The bilinear quark currents in the original 3d Euclidean theory takes the following 
form, for the vector current V"", 

t>3 = -T^e-*('^"-'^'")"V'mV^n (2.7) 

m,n 

= eTEe'*^""""'"^"?„^7'^n (2.8) 

m,n 

with j = 1,2 and ^ = z^^"^). 

We now state some key results for the 't Hooft model that we will use in the following. 
For a review and references to original papers, see e.g. [^. Due to the linear Coulomb 
potential V{r) = or, with the string tension a = g'^Nc/ir, there are no free quark states, 
but only color neutral hadrons, and at energies much larger than the quark masses, the 
bound-state spectrum is approximately given by 

Ml = fila = Ti^ak A; > 1 . (2.9) 



In real QCD, the Greens functions for the vector and axial vector currents are of 
special importance, since they couple to external electro-magnetic and weak probes. In a 
model problem one is of course free to consider any arbitrary Greens function, but we shall 
concentrate on the vector current two point function, i.e. the color singlet polarization 
tensor, li^yiu^q), that at T = determines the (would be) e+e~ production rate, at 
T 7^ and u; = gives the screening of external electric and magnetic fields, and for 
oj > \q \ contains information about the finite T quasi particle spectrum. At T = this 
function has been extensively studied, and, in the case of one quark flavor and to leading 
order in l/N,,, the following exact expression has been derived [0, [2^, |21|, ^ 



where 



f7f)= f'dxM^) , (2.11) 
Jo 

and 0„(x) is the meson wave function and M| its invariant mass. Note that the mesons 
have zero- width to leading order in 1/Nc. 

A similar analysis for the scalar current two point function in the leading approx- 
imation yields 

n.„^).-(^]i:^M^ . (2.12) 

where fi, fin is the four momentum of the current and the bound state mass, in units of 
g'^Nc/ 71, with the coupling of the the bound state to the respective source given by 

1 /"^ , / \/^ \FTh 



7a, 76 are the quark masses, also in units of g^Nc/r^. However, due to parity, 

ma f — 0„(x) = i-ir+'m, C <Pn{x) , (2.14) 

Jo X Jo L — X 



SO the sum over the bound states in eq. (|2.12| ) is effectively restricted to even n. One can 
also show that, asymptotically, 

^?f^-^ • (2.15) 



Similar relations can easily be derived for the pseudoscalar current. 

For nonstatic correlation functions the situation is more complicated. First, we must 
extend the above treatment to include bound states of quarks with different masses. 
To do this we have generalized the formulas ( p.lO|) and (|2.12| ) to the case of unequal 
masses, i.e. to couplings of the form ipnlf^i^m where the field ipn has mass M„. The 
resulting expression for the polarization tensor is identical to ( p.lO| ), but with the wave 
functions being solutions to the unequal mass 't Hooft equation; the same holds for the 
scalar correlator ( 2.12|) . Although the result is simple and intuitive, the derivation again 
involves several nontrivial cancellations. There is however a more subtle issue that regards 
the nonstatic gluon modes. As seen from ( |2.3| ), they manifest themselves in the 2d theory 
as heavy vector fields in the adjoint representation. Since there is no physical light gluon 
in two dimensions, we expect these heavy vector particles to be confined, with a string 
tension aAd = (Caci/Cf)'^, where C^d and Cp are the Casimir operators in the adjoint 
and fundamental representations respectively. For a general discussion of 2d Yang-Mills 
theory with adjoint matter fields, see e.g. [l^. In the large Nc limit we expect 

asymptotically straight trajectories of stable bound states of the type ipj^_^Ami....Am^.4'n2j 
and there is no reason to expect that they should not contribute to correlation functions. 
In fact we would expect a general formula of the form 

where and are masses and couplings for these trajectories. An important point, 
to be discussed below, is that Lorentz invariance requires the vector particles to couple 
either in pairs or via derivative couplings. 



3. Static and dynamical quark susceptibilities 

In the large Nc limit, all contributions to color singlet correlation functions will be of 
the form given in (|2.10 ) and ( |2.12 ). Using the results from the 't Hooft model, we can 



calculate the contribution from the qq trajectory to the finite temperature correlator of 
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vector currents, 



G^"'{nN,q) = dr f (fre-'^^^''^'^'^{V>'{T,r)V''{0,0)) , (3.1^ 







where the expectation value is with respect to the action S^e- We have set f2jv = '2'n'NT 
and used r = {x,y). Using (|2.6| ) and ( p.7|) we can relate G^^ to the 2d scalar correlator, 



and by combining this with (|2.12| ) we get 



, (3.2) 

where M^{m,n) is the mass and (7™'", given by (p.l3| ), the form factor of the k:th pole 
on the trajectory describing mesons consisting of quarks with masses {2n + l)7cT and 
(2m + 1)itT. Similar expressions can be given for correlators of other currents. Note that 
since we have rotated to a 2d Minkowski space, = Iq—Qi- After having evaluated ( p.2| ) 
using 2d QCD methods, we can obtain the corresponding real time finite T correlation 
function by the analytic continuation iQ^ = 00 + ie. 

As shown in appendix A, perturbation theory to one loop gives the following leading 
behavior at large T and \q \ = q u, 

G^'{uj,q) = -^(T\n2-^)+0{qyT) for T > g > ^ . (3.3) 



TT V 16. 

This expression is also leading in Nc, so we can hope to reproduce it using the 't Hooft 
model, and we shall now demonstrate that it can be calculated from (|3.2|) , i.e. by neglect- 
ing both the contributions from the bound states containing heavy scalars, and the effect 
of the light Higgs field. The derivation is a straightforward extension of the calculation 
in which we now briefiy recall. Since each term in the sum (|3.2| ) is ~ 0(1) in the 
temperature, a result ~ T must emanate from the infinite range of the summation. Thus 



we can use the asymptotic expressions ( p.9|) and (|2.15|) for bound state masses and scalar 
form factors, and furthermore replace the discrete sum with an integral of a continuous 
bound state spectrum, 

E-i/ , JM^ > (3-4) 



k 



TT^CT J(mi+m2) 



with mi and m2 the two different quark masses. The resulting expression for the relevant 
correlation function (after the rotation —>■ — g^) reads, 

where we introduced an ultraviolet cutoff A. The static limit corresponds to = 0, 
where ml = m\ = u"^. Performing the M^-integration and rewriting the sum over Un as 
a countour integral we get for the T-dependent piece, 

Gg(^,g) = ^n,m\n f^^) , (3.6) 

where np is the Fermi distribution function. The logarithmic function has a cut along 
the real axis; in the positive real half-plane the cut starts at q/2. Following |18[ we can 
deform and close the contour to get 



G|:(c. = 0,g) = -^(Tln2-|) (3.7) 



^33 



and comparing this with (|3.3| ) we see that the leading terms agree, but the subleading 
ones deviate, although they have the same fc, T dependence and the same sign. 

We now turn to the non-static case. In appendix A we derive the following perturbative 
result, 



Gf,{uA) = ^^ln2[-^)+0(VM forr»cu»g . (3.^ 



By comparing (|3.^ ) with (|3.3| ) we see that the two limits g = , u; — > and u = , q ^ Q 
do not commute, just as in 3+1 dimensions. 

The limit (|3.8| ) is a bit harder to reproduce with the 2-dimensional technique than 
the static one, since mi 7^ 1712 which implies a more complicated analytic structure in the 
2;-plane. Performing the M^-integration we get 



^i(-,g) = -^/— tanh^ 



g2+ V-(^-^)2 



(3.9) 



where we analytically continued back to real time by zfi^v = w + ie, and where the 
integration contour C and the square root cuts in the z-plane are shown in fig. 2(a) in 
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appendix B. After some exercise in complex analysis, which is outlined in appendix B, we 
get 

G33(^, q) = ajNc —npiPx) fin (q^ + n% - Ax^ - Aixnj,) 

- ln(g2 + ^]^-4x2 + 4ixfiiv)] ■ (3.10) 



Taking g = in ( |3.10|) the logarithms cancel for = 0, so we conclude that in the 



high-temperature limit, 



GfJuj,6) ^uj + 0{l/uj^) forcu<T . (3.11) 



Note that the non-analytic properties of the perturbative results are also manifest in ( |3.9|) , 
since if we take the limits in reverse order we get the logarithmic cut leading to ( |3.7|) , 
instead of (|3.11|) . We have thus reproduced the nonleading behavior for the perturbative 
result (|3.8|) , in the limit q = 0, but from the expression (|3.10 ) it is also clear that for 



q the leading term will not pick up any linear T-dependent piece, but only 0{q^ /uj). 

At the end of the previous section, we pointed out that there is no compelhng reason to 
neglect the contribution from trajectories containing heavy vector particles. We now want 



to argue that the leading term in ( |3.8| ) is a likely sign of such a trajectory. The simplest 
possible contribution would come from bound states of the type Tp^^An'ipm2 illustrated 
in fig. 1. Since A„ carries a vector index the coupling must be derivative, i.e. of the 
form gqAqq^i where the coupling constant QqAq has dimension 1/mass. Without solving 
the resulting 3-body problem we can not calculate this coupling, but we can nevertheless 
predict that the contribution to the correlation function from the trajectory with mi = 1712 
and m^v = f^AT = 2ttNT will be 

Gg,(^,gl=^?|A,^^Gg(0,g) = -^|^/— ^ln2 . (3.12) 

Although we have not determined the coupling, the presence of the factor Tin 2 in this 
formula is very suggestive when comparing to the perturbative result ( pj.S] ). It is very 
natural to expect it to have the same origin, i.e. from an infinite summation over states 
on a linear trajectory. For (|3.12|) to be consistent with ( |3.8D we must have g-gAq ~ ^/^n, 
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i.e. the scale is set by the mass of the scalar particle. We have not been able to prove 
this. 



Fig. 1. Contribution to the bound states containing a heavy vector particle. 

To summarize this section, we have found a rather pleasing picture in terms of the 
interplay (high temperature duality) between perturbation theory and the dimensionally 
reduced theory at high temperature. In the static sector only the lightest gluon modes 
need to be retained in order to correctly reproduce the screening length, whereas the 
non-static sector for consistency requires bound states containing heavy vector particles. 



4. DeTar correlators at small q. 

We argued above that the formula ( p.2|) can be used to calculate the non-static DeTar 
correlator in the limit q = 0. However, to be able to make the final analytical continuation 
Qn = —i^ we need expressions for both bound state masses and form factors that are 
analytic in the masses, i.e. in fijy- Clearly some approximation is needed, and we have 
tried to use the semiclassical methods developed in Even this turns out to be too 
difficult, since the WKB expressions for the wave functions in general are very complicated, 
and do not allow us to derive explicit analytic expressions for the form factors. However, 
the nonrelativistic system of two particles with masses mi and 1712 moving in a linear 
potential V{r) = ar is relatively easy to handle. The bound states of this system are 
given by, 

2/3 

Ml{m,,m2) =M^ + 2\kMj^-^ + 0(mj/') (4.1) 

for a <^ nil , ^2, with M = mi +m2, m = mim2/M and A„ the n*^ root of the Airy func- 
tion Ai'(x). We can also calculate the nonrelativistic form factors which are proportional 
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to the wave function at the origin, 

{gT^-^'f - - {maf" (4.2) 

for even parity states. For odd parity states, the wave function vanishes at the origin. 
Although in principle straightforward, we shall not elaborate on the evaluation of the 
normalization in ( [4.2|) . We are interested in the result after the analytic continuation 
f^Tv = —ioJ^ when the masses become, 

mi = |2n + l|7rT (4.3) 
m2 = + {2n + l)vrT| , 

where |x| = (x^)5 which implies that Re(mi + m2) = \2n + l|27rT. In particular, we get 
for the lowest masses, 

^^-2^^±^^+(2^(^±(^T)^j ' ^'-'^ 

where the upper and lower sign refers to n = and n = —1 respectively. 

The question now is whether we can use this information to calculate the correlator 
( p.l|) , at least in some specific kinematic region. This is not a simple problem. First notice 
that the correlation function cannot be calculated directly for fixed u since the sum over 
the radial states does not converge fast enough for the nonrelativistic approximation to be 
useful. Also, the analytic structure in u obtained from ( p.lO|) and ( |3.12|) is too complicated 



to allow us to evaluate the real time correlator by a fourier transformation. The way to 
proceed is instead to consider the following Fourier transform, 

G"^^(fijv,x) = 1 1 ^e*^^-'G^^(n^,0,g2) (4.5) 
and calculate the contribution from the gg-trajectory only, to get 

^2 AT oo 

Gil{^N,x) = -z^-^ E j:i9k^'''yMk{n + N,n)e-^^^^-+''^-> . (4.6) 

n=-oo k 



In the limit xT ^ 1 and <^ T, the expression (i^) is saturated with a few nonrelativis- 



tic bound states. Note that in the static limit = we retain the static DeTar correlators 
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derived in 0, and the main result of this section is that the x— dependence of the dynami- 
cal correlators are obtained by a straightforward extension of the earlier results. Also note 
that large x corresponds to small q. Following the arguments of the previous sections, we 
expect trajectories containing heavy bosons to give only small contributions. 

We shall now make the conditions that allow us to saturate ( |4.6|) with only a few 
bound states more precise. As already stated, depends linearly on |n|T, so for xT ^ 1 
equation ( [4.6| ) will be dominated by the two lowest masses, corresponding to ([4.4|) . To 
suppress higher radial excitations we must require 

x»Q' , (4.7) 

which combined with the weak coupling condition T ^ (7 ~ y/a, implies 

xT>|^^y»>l , (4.8) 
which means that we are in a region where the correlation function is very small. 



5. Discussion and Conclusions 

We have shown how to relate real time finite T correlation functions in 3d QCD to 
quantities that can, at least in a certain limit, be calculated analytically in a generalized 
't Hooft model with an infinite number of heavy quarks and scalars. Several assumptions 
and approximations went into our calculation and we want to comment upon their validity. 

First, we neglected the Higgs field. This is a questionable step, and might well be 
improved upon. The effect of the Higgs field can be estimated in two limits. For the 



stringlike excitations far up on the "Regge" trajectories we can use the arguments of [|T7 
to conclude that the only effect of the Higgs field would be to renormalize the string 
tension and add a non- leading 1/R term to the potential. This will clearly not change 
the leading behavior of the spectrum which was the only thing needed to retain the 
perturbative result for the susceptibility. For the nonrelativistic bound states we can 
estimate the effect on the masses from Higgs exchange to lowest order in perturbation 
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theory, using an effective four fermion potential of tlie kind discussed in . However, as 



long as there are no numerical simulations to compare with it is hardly worth elaborating 
on corrections of this type. 

Secondly we want to emphasize a point already made in section 3. Starting from 
the high temperature limit, we showed that consistency between the perturbative high 
temperature theory [i.e. 3d QCD) and the generalized 't Hooft model requires that the 
spectrum of the latter contains trajectories of bound states with heavy vector constituents, 
and that their decay strength decreases as l/u. These effects are dominant for time-like 
separations, whereas e.g. the screening length can be calculated by only taking into 
account the gg-trajectories. 

The limited results of this paper may be compared with the very detailed ones previ- 
ously obtained in the Schwinger model and the large limit of the 2d Gross-Neveu 
model where the analysis is simpler. The Schwinger model is exactly solvable in terms 
of an elementary scalar field, and the 2d Gross-Neveu model is tractable in the ap- 
proximation. In contrast, 3d QCD studied here is rather similar to real QCD, and the 
difficulties are accordingly increased - in fact we find it encouraging that nonperturbative 
statements about the asymptotics of the Detar correlator, admittedly in some specific 
kinematical domain, can at all be made using analytical methods. 

Finally we want to comment upon to which extent the results of this 3d model can 
be generalized to real, (3+l)d QCD at finite temperature. Starting from a high T phase, 
the correlators in 4d QCD show a similar non-analytical behavior as encountered in 3d; 
in the static case we then have G|g ~ T^, and in the nonstatic, ~ T^(g^/co'^). Since 
it is also known that QCD in 3d is confining, it is from this point of view not unlikely 
that the dynamical correlators in the generalized 3d theory have to be accompanied by 
a summation over trajectories of bound states consisting of heavy vector fields as well, 
with a derivative coupling and a coupling constant g-gAq ~ l/t^, just as in the 2d case. 
However, since 3d QCD is not solvable even in the large limit, this conjecture is not easy 
to verify. For the dynamical DeTar correlators, we expect the asymptotic x-behavior to be 
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generic, i.e. that the sum over the bound states is saturated by only a few nonrelativistic 
states and that the leading u dependence of the correlation functions comes from the 
uj dependence of the screening masses just like in ( [4 .41) and ( [4.6|) . Again there will be 
corrections from the mass-dependence of the binding energy. 

For non- asymptotic x, nothing can be said for sure, since many bound states will con- 
tribute, and the nonrelativistic approximation can no longer be trusted. One should note, 
however, that to the extent that the constituent masses enter the total mass additively, 
the leading u dependence in ( [4.6|) could still be correct, since the factor e**^^ is common to 
all the terms in the sum. If valid, this argument, which is the same in both three and four 
dimensions, predicts a very simple u dependence for the DeTar correlators, and might 
even be related to dynamical quark-gluon plasma parameters. Again we stress that these 
comments on the non-asymptotic behaviour are highly speculative. 



Appendix A. The polarization tensor at high T 



In this appendix we give the derivation of the asymptotic results, equations ( |3.3| ) and 
(p. 81). The polarization tensor G{Qn,(i) = {Ncg'^)^^G^^{QN,(l) is given by 



. (Pp (I 1 
^2^ / T^r^Tr I as-cTs 



(fp fpl + p'^+Poqo+p-q 



(A.l) 



(27r)' V P^{p + q)^ J ' 
where po = (2n + l)7cTi and go = ^NrcTi. Evaluating the frequency summation with 
the help of contour integration, omitting the vacuum piece and performing the analytic 
continuation iflj^i = u + ie, we find for the imaginary part 



ImG(ti;, q 



-e{uj - q) 



dp 
2^ 



4p2g2 _ ^^2 _ q2 _ 2pujy 



{npiPip - uj)) - npiPp)) 



dp 



4p2 _ _|_ ^2 _ g,2 

4p2g2 _ ^^2 _ g2 _ 2puj)'^ 



(A.2) 



where 6 is the step function. The integration limits are g+ = {q+uj)/2, uj±. = {uj±q)/2 and 
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9 = 1^1) P = \p\- From this expression it is clear that the imaginary part is proportional 
to u when q ^ u and to g for a; q. 

The real part of the polarization tensor is 

T f<^ dp , r, \ i + 4pco' + tu^ - 



1 r°° (IT) 

-ReG(cu,g) = -ln2-0{uj~q) —UpiPp) 

TT JO 2tI 



'uj^ — q^ + 2pujY — ^p^q^ 



. dp V-4pa; + w2-g2 

+e{iu-q) / —nFiPp) 



k+ 27T (^^2 - q2 _ 2pijjf - 4p2g2 

4p2 _ _|_ _ g,2 

(u;2 — q^ — 2pLoY — ^p^q^ ^ 

4p2 _^ _^ ^2 _ ^2 

— g2 + 2pujy — Ap'^q^ ^ 



Jo ZTT \ 



•^0 27r W(cu2 - g2 _ 2pcu)2 - 4pV 



Expanding in powers of q/uj, we find 

27r lcu2 



ReG'(^,g)^^f4Vn2 (A.4) 



in the limit T ^ uj ^ q. 

In the other limit, T ^ q ^ uj, the expansion of ( [A .31) gives 



ReG(u;,g') = In 2 + -^ . (A.5) 

vr lb 

Equation ([A.5|) gives the screening mass 

m^; = -(732G(a;,g = 0) = (?32-ln2 , (A.6) 

IT 

a result that also follows from the pressure of the free quark-antiquark gas, with a non- 
vanishing chemical potential. 



Appendix B. The polarization tensor from 2d, q = 

Here we show that the polarization tensor calculated from the 2d theory has no leading 
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T-behavior for q = 0. Starting from the expression (|3.5| ) and performing the M 



integration, we can rewrite the Matsubara summation as 
1 



+ {y/^z^ + \/-(z-in 



N 



(B.l) 



where the contour C is as shown in fig. 2(a). We have two different cuts for the square 
roots, shown as wiggly hues (for definiteness we have chosen Qn > 0); the cut for z = i^N 
is always at an even multiple of inT and hence it will not interfere with the poles of 



tanh(/?2/2). These cuts correspond to having the cut for \j f{z) along the negative real 
axis in the complex /-plane. 







i 









(a) (b) 
Fig. 2. The integration contour for G'^^, before (a) and after (b) deformation. 

We can now deform the contour and arrive at the contour shown in fig. 2(b), where 
we have neglected T-independent pieces; this follows from the identity 

1 



1 , PA , 1 

- tanh — = ±- =1= , „ 

2 12/ 2 e±^^ + 1 



(B.2) 



The cut at z = then gives 



and by using ( |B.2| ), equation ( [B.3| ) becomes for the temperature dependent piece 



— / —npiPx) [in (-4x^ - AixVLn + + q 
Itx Jo 2711 L ^ 

- In (^-4x^ + AixVlN + fi^ + g ^ 



(B.4) 



In the same way, we find for tlie second cut at z = ifipf, 



1 r°° dx 
2^ 



Jo 27ii L ^ 

- In (^-Ax'^ + Aixn^ + n% + q^)] (B.5) 



by using 



tanh tanh (B.6) 

and so equation ( |B.1|) becomes 

G(QN,q) = - npiPx) In (Q'i + q'^ - ix^ - Aix^N] - 

Tc Jo 27ri L ^ ^ 

- In (n% + q^ - Ax^ + AixVlN)] ■ (B.7) 



Taking g = we see tliat by putting Qn = ^ tlie integral vanishes. Since it also has to 
vanish for /5 — > cxd we conclude that, after analytic continuation, 

G33(^,g = 0)~a7exp(-^) , (B.8) 

which is a subleading term and goes like ~ u for high temperatures. 
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